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ABSTRACT 

We update our prior work on the case B collisional-recombination spectrum of He I to 
incorporate ab initio photoionisation cross-sections. This large set of accurate, self-consistent 
cross-sections represents a significant improvement in He I emissivity calculations because 
it largely obviates the piecemeal nature that has marked all modern works. A second, more 
recent set of ab initio cross-sections is also available, but we show that those are less consistent 
with bound-bound transition probabilities than our adopted set. We compare our new effective 
recombination coefficients with our prior work and our new emissivities with those by other 
researchers, and we conclude with brief remarks on the effects of the present work on the 
He I error budget. Our calculations cover temperatures 5000 ^ T e ^ 25000 K and densities 
10 1 < n e < 10 14 cm . Full results are available online. 
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1 INTRODUCTION 

The high precision required to make a significant measurement 
of the primordial helium abundance presents challenges that are 
unique in nebular astrophysics. The recombination rate coefficients, 
needed to convert emission line intensities into ionic abundance ra- 
tios, must have an accuracy better than the precision expected for 
the derived abundance. Usually we hope to measure abundances of 
most elements to an accuracy of 20 - 30%. The primordial helium 
problem requires recombination rates accurate to better than a per- 
cent, presenting unprecedented challenges to the atomic physics of 
the line formation. See Brocklehurst (1972) and references therein 
for early seminal work in the field of theoretical He I emissivities. 

We need four types of atomic data to calculate the case B re- 
combination spectrum: level energies, transition probabilities, pho- 
toionization cross-sections, and collision rates. Porter et al (2009) 
summarize the contributors to the error budget. Energy uncertain- 
ties have always been negligible compared to other sources, and 
accurate transition probabilities have been available since Kono 
& Hattori (1984). Drake (1996) and Drake & Morton (2007) im- 
proved upon the latter still further. 

Photoionisation cross-sections and collision rates represent 
the greatest sources of uncertainties in current standard calculations 
of He I emissivities (Benjamin, Skillman, & Smits 1999, hereafter 
BSS99; Porter et al. 2005, 2007). For the low-density, extragalactic 
observations used in primordial helium analyses (Peimbert, Lurid- 
iana, & Peimbert 2007; Izotov, Thuan, & Stasinska 2007; Aver, 
Olive, & Skillman 2010), photoionisation cross-sections (and, by 
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extension, recombination coefficients) are the greatest source of un- 
certainty (Porter et al. 2009). See Ferland et al. (2010) for a recent 
review of the errors in He I emissivities and in primordial helium 
abundances. 

In this Letter we update our earlier calculations of case B, He I 
emissivities (Bauman et al. 2005; Porter et al. 2005, 2007) to in- 
clude a large set of self-consistent, ab initio photoionisation cross- 
sections. We compare our present results with our previous work, 
with those by BSS99, and with those by Almog & Netzer (1989). 
We present here a subset of our present results that is most applica- 
ble to primordial helium research. A much larger set is available in 
the electronic edition. 



2 ATOMIC DATA 

The first large set (L ^ 4 and n ^ 25) of modern, self-consistent 
photoionisation cross-sections was calculated by Hummer & 
Storey (1998, hereafter HS98). However, the full set has remained 
unpublished. Ercolano & Storey (2006) used the full set to produce 
accurate continuous emission spectra. The current work represents 
the first use in case B line emissivity calculations. 

Another large set of photoionisation cross-sections was pre- 
sented recently by Nahar (2010). Because highly-accurate bound- 
bound absorption oscillator strengths obtained with Hylleras-type 
wave functions are available (Drake 1996), we can readily apply a 
simple extrapolation method (Burgess & Seaton 1960) to judge the 
accuracy of free-bound differential oscillator strengths at the ion- 
ization threshold. HS98 applied this method to several low-lying 
levels of He I and tabulated results in their Table 1. In Figure Q] we 
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Figure 1. Percentage differences for threshold photoionisation cross- 
sections of ns 'S states relative to "extrapolated" Drake (1996) results. Nahar 
(2010, blue circles), HS98 (red crosses). 



Figure 2. Ratio of present to PFM07 effective recombination coefficients in 
the low-density limit at 10, 000K. The three levels with the largest changes 
are 3d [ D, 3d 3 D, and 5i'S. See text. The differences will generally be 
smaller at finite densities. 



compare those results with the ab initio calculations of ns S lev- 
els from both Nahar and HS98. Nahar's values disagree with the 
extrapolated cross-sections by as much as 4%. The HS98 cross- 
sections are clearly more consistent with the Drake results than are 
the Nahar (2010) values. The target wave functions used in HS98 
were derived from those of Calvert & Davison (1971) and were 
designed to fully represent the dipole and quadrupole polarizabil- 
ity of the He + ground state and to include short-range correlation 
optimised on the energies of the ns l S states, which were known 
to be particularly difficult to represent. This target provides very 
accurate results at energies near threshold, as required for the tem- 
perature range and applications discussed here but is unsuitable at 
higher energies. The target used by Nahar (2010) is extensive and 
implicitly treats polarization and correlation effects but it is not op- 
timised for near threshold calculation or for the ns S states. It is, 
however, applicable over a much larger energy range than that used 
in HS98. 

We use the unpublished cross-sections described by HS98. For 
n ^ 4 and L ^ 2, we rescale to agree at threshold to those extrap- 
olated from the Drake oscillator strengths and listed in table 1 of 
HS98. We use the extrapolation procedure described by HS98 to 
obtain values for n = 5 and L ^ 2, but these normalisations have 
only minor effects on the emissivities of the most important lines. 

We follow HS98 and assume that recombination coefficients 
are fixed relative to hydrogen for n > 25 and L ^ 2. All other aspects 
of our calculations are as described in Bauman et al. (2005). Briefly, 
this entails rcLS-resolved levels for n ^ 100, a single "collapsed" 
level at n = 101, and collisions induced by electrons, protons, and 
He + . The calculations are performed with a development version 
of the spectral simulation code CLOUDY (last described by Ferland 
etal. 1998).Q 

Following Switzer & Hirata (2008), we also have added radia- 
tive bound-bound electric quadrupole transitions (Cann & Thakkar 
2002). However, these transitions have a negligible effect on the 
present results; we mention the change for completeness and to 
note that these transitions may be important for high-Z He-like ions 
(as in Porter & Ferland 2007). 

1 The /-resolved code described by Bauman et al. (2005) has also been 
updated with the HS98 cross-sections. That code can be retrieved with a 
Subversion client from https://svn.nublado.org/bauman/source. 



Finally, we note that we have compared our existing imple- 
mentation of heavy particle angular momentum changing collisions 
to the new results by Vrinceanu, Onofrio, and Sadeghpour (2012) 
and found excellent agreement. 



3 RESULTS AND DISCUSSION 

We compare effective recombination coefficients in our new treat- 
ment with those of Porter et al. (2007) in Figure|2] For the majority 
of levels the change is < 1%. There are two clear exceptions, corre- 
sponding to 3d l D and 3d^D. The cause is a programming error in 
our earlier renormalisation of Peach (1967) cross-sections. There 
are only two levels affected. Unfortunately, these are the upper lev- 
els of two of the most important lines in primordial helium abun- 
dance works: AA5876 and 6678. The new results yield stronger 
emissivities for both lines. The next largest differences in Figure [2] 
correspond to levels with only weak optical lines. 

We compare our new emissivities with BSS99 for several 
strong lines in Figure[3] This figure is directly comparable to figure 
5 of Aver, Olive, & Skillman (2010). Consistent with the discussion 
above, j^76 an d 76678 are now m better agreement with BSS99, 
though important differences clearly remain. 

Our calculations cover 5000 < T 25000 K (in 1000 K steps) 
and 10 1 «S n e sC 10 14 cm (in 1 dex steps). In addition to the lines 
published in our prior work, we have also included several weak 
infrared lines. Table [T] lists the wavelength and upper and lower 
level designations of all reported lines. Table [2] contains a small 
subset of the full results that is most applicable to primordial helium 
abundance calculations. 

Almog & Netzer (1989, hereafter AN89) also presented He I 
emissivities for electron densities up to 10 14 cm~ 3 . In Figure [4] 
we compare present emissivities in AA4471, 5876 and 7065 to the 
AN89 results (with = 0). The calculations generally agree for 
n e < 10 8 cm~ 3 (where collisions do not dominate over radiative 
processes) and for n e ~ 10 14 cm~ 3 (where collisions are so domi- 
nant that levels are approaching local thermodynamic equilibrium). 
Critical densities for levels with n ~ 3 or 4 fall between these re- 
gions. Collisions are dominated by excitations from the metastable 
2s^S, which is treated in both the present work and AN89. Our 
work clearly has enhanced collisional excitation relative to AN89, 



© RAS, MNRAS 000,[T]E] 



Improved He I Emissivities 3 



1.05 | r 

3889 

1 04 - 4026 
1 U4 4471 

5876 

1.03 - 6678 

7065 

I 1-02 = 



1.01 



0.99 



0.98 



0.97 



Table 1. Wavelengths and upper and lower levels of reported lines. 
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Figure 3. Ratio of BSS99 and present emissivities for several strong lines 
as a function of temperature with n e = 100 cm~ 3 . This figure should be 
compared with figure 5 of Aver, Olive, & Skillman (2010). The principal 
effect is that for both A A5876 and 6678, the present results are greater than 
the earlier Porter et al. works and are roughly the average of the BSS99 re- 
sults and our earlier work. Our A A3889 and 7065 results have also changed 
considerably, and (within the plotted temperature range) the ratios no longer 
cross the "zero deviation line" discussed by Aver, Olive, & Skillman (2010). 
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Figure 4. Comparison of present emissivity calculations to those by AN89 
as a function of density at 10.000K. 
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and differences between the two calculations are strongly corre- 
lated with the relative collisional contributions given in table 5 of 
Porter et al. (2007). Note, however, that those collisional contribu- 
tions should not be added to the present tabulations; they are al- 
ready included. 

A recalculation of the Monte Carlo error analyses performed 
in Porter et al. (2009) is beyond the scope of this paper. We expect 
that our analysis of the relatively high-density 'Galactic' model 
would be largely unchanged, as collisional uncertainties should 
dominate both before and after the present work. In the low-density 
'Extragalactic' model, however, our 'optimistic' uncertainties now 
seem somewhat more realistic. Total uncertainties are largely due 
to uncertainties in threshold photoionisation cross-sections. One 
measure of those uncertainties (used by Porter et al. 2009) is the dif- 
ference (< 0.7%) between the extrapolated threshold cross-sections 
and the HS98 ab initio results. Goodness of the fits used to extrap- 
olate threshold cross-sections, including sensitivity of the fits to the 



number of terms used in the fitting procedure, is another measure, 
and we obtain < 0.2% for all levels with n ^ 5. 
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Table 2. Emissivities of several He I lines at conditions important for primordial abundance analyses. This table is a small subset of the full results. Values are 
Allj '/« c mhc+ in units 10~ 25 erg cm 3 s _1 . 
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